Impurity and edge roughness scattering in armchair graphene 
nanoribbons: Boltzmann approach 
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Abstract -The conductivity of armchair graphene nanoribbons in the presence of short-range 
impurities and edge roughness is studied theoretically using the Boltzmann transport equation 
for quasi-one-dimensional systems. As the number of occupied subbands increases, the conduc- 
tivity due to short-range impurities converges towards the two-dimensional case. Calculations of 
the magnetoconductivity confirm the edge-roughness-induced dips at cyclotron radii close to the 
ribbon width suggested by the recent quantum simulations. 



Introduction. — The theoretical descriptions of elec- 
tron transport in graphene sheets with various scattering 
sources, such as charged impurities in substrates, micro- 
scopic corrugations or short-range resonant scatters, have 
been frequently based on the Boltzmann approach. Such 
studies, in particular the electron density dependence of 
the conductivity, are of fundamental interest since they 
help identifying the dominant scattering sources. [TJ[2] 
It has been verified systematically that the Boltzmann 
approach works quite well to describe the transport in 
broad parameter ranges for both single and double layer 
graphene. [3] As the width of the graphene strips is de- 
creased, graphene nanoribbons (GNRs) are fromed, where 
size-dependent effects, for example the inhomogeneous 
electron density [4] or the edge roughness, become rele- 
vant for the transport properties. In this case, the trans- 
port has so far been described within the framework of 
the Landauer-Biittiker model with the aid of Green's func- 
tion techniques. [51|B] It is well established that in GNRs, 
edge disorder can contribute significantly to the scatter- 
ing [THlOj which in wide structures is governed by a com- 
bination of scattering at charged impurities and resonant 
scattering at short-range defects. Edge disorder has been 
suggested as the source of the transport gap in narrow 
GNRs around the charge neutrality point. [TU1 - [T5] Further- 
more, a typical size effect, the so-called edge-roughncss- 
induced magnetoconductance dip (ERID) in GNRs has 
been studied by numerical quantum simulations, which is 
interpreted as a magnetic-ficld-cnhanccd diffusive scatter- 



ing when the electron trajectory grazes at the edges. [14] 
On the other hand, the Boltzmann approach has been ap- 
plied to treat a variety of scattering sources in conven- 
tional quasi-one-dimensional (QID) systems, for example 
quantum wires. |15fll7| . However, only a few aspects of 
transport in GNRs have so far been studied within the 
Boltzmann model. [TSlfTO] 

In the present paper, we apply the linear Boltzmann 
equation to armchair GNRs and determine its transport 
properties in the presence of (5-type short-range impurities 
and edge roughness. The magnetoconductivity in wide 
GNRs with rough edge roughness is studied as well. 



Model and theory. 

Hamiltonian 



We start with the Dirac 



(1) 



with Fermi velocity vp ~ lO^m/s and Pauli matrices 
ax.y and acting on the A/B sublattice and K/K' val- 
ley spaces, respectively. The energy spectrum of GNRs 
depends on the nature of their edges, namely zigzag or 
armchair. Within the present work, we restrict ourselves 
to metallic armchair GNRs (AGNRs). For this system, 
the boundary conditions imposed on the wave function, 
namely ^'^(a; = 0) = "isix = 0) = ^a{x ^ W) ^ 
"^six = W) = 0, give rise to the allowed transverse wave 
vectors as 
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where fc^ is Fermi wave vector in the Jith subband. The 



with a ~ 0.246nm being the lattice constant of graphene 
W is the width of AGNRs. The integer n is of the order transition matrix element Tnn' is defined as 
of W/a for the energetically lowest modes. Throughout 
this text, we denote the energy e normalized to Tivp as y 
= e/{hvp) with — + . The normalized wave " 
function for the ?7th subband reads ppll^ 
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with the summation running over the mode index fi. The 
Boltzmann conductivity for GNRs is then given by 



which is a mixture of two Dirac points K = (47r/(3a), 0) = 
(K, 0) and K' = (-47r/(3a), 0) = {-K, 0). Here we choose 
the wave vectors in the x-direction to be quantized and the 
transport is oriented along y-direction. L is the length of 
the system. 

To describe the transport properties of GNRs, we adopt For nonzero temperature, the conductivity is obtained 
the linearized Boltzmann equation describing the general from 



h nEp W 
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QID system 
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dk 
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We proceed by describing how the scattering potentials 
have been implemented in this formalism. First, we con- 
sider (5-type impurities in the form of 



Ni 



where Ey is the applied electric field along the transport 
direction, fnk is the distribution function of a state with 
wave vector k and energy e„fe in the nth subband, and the 
superscript "0" denotes the equilibrium distribution. Ac- 
cording to Fermi's Golden Rule, the scattering probability ^hcre 7 and Nj are the strength and the number of im- 



U = lJ2^^^~ - yj) 

i=i 
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due to the perturbation potential is given by 

W„'fc'„fc = —\{n\k'\U\n,k)\'^S{en'k' - f-nk) 



purities, respectively. Thus the matrix element squared of 
the perturbation is evaluated as 



(5) 



\{n'k'\U\nk)\'' = 
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Using the relaxation time approximation, the noncquilib- 
rium distribution function can be written as 

/«fc(e«fc) = fnki^nk) - e-B^w„(e„fc)T„(e„fc)5(e„fc - Ep) (6) 

with Fermi energy Ep and the relaxation time t„ for the 
state in the nth subband. The velocity for the nth sub- 
band Vn = {l/h)denk/dk = vpk/^y (kn)'^ + fc2 for the lin- 
ear spectrum of graphene. Inserting Eq. ([6]) into Eq. (jlj, 
the Boltzmann equation at zero temperature can be writ- 
ten as 
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with Hi ~ Nj/WL. The transition matrix elements are 
finally given by 
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(7) Using Eq. ^ and Eq. ^ we obtain for the conductiv- 
ity of GNRs the expression 



Multiplying both sides of Eq. ^ hy k and summing over 
fc, we obtain after some algebra ^ 
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For a large number of subbands, i.e., A/" 3> 1, we have 
kfj,,kn <C Ep and kp,kp ^ Ep. In this approximation, 
Eq. ^T5\i converges to the well-known result for the ex- 
tended case 



8e^ 
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which is independent of the carrier concentration. 

As a second scattering mechanism, we study the effects 
of edge roughness in the absence of magnetic fields. The 
edge roughness is parametrized by A{y)dV{x) /dx, an ex- 
pression which has been applied before to model rough 
semiconductor quantum wires and interfaces. j22ll23j V{x) 
is the one-dimensional confinement potential which can be 
modeled by a finite mass term in the Dirac Hamiltonian. 
[24] A{y) is a function describing the potential fluctuation 
of GNRs and characterized by (A(y)) = and the auto- 
covariance function {A{y)A{y')) = exp[— (y — j/')^/A^] 
with A being the correlation length. Furthermore, (• • •) 
denotes position averaging. 

To evaluate the perturbation matrix element, we define 
the function S related to the a;-componcnts of the wave 
functions as 
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where 4>n{x) denotes one of the components of wave func- 
tion in Eq. ([3]). For the hard- wall confinement potential 
present in GNRs, this function can be expressed as 
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It is noteworthy that a linear form for the matrix ele- 
ments of the edge roughness perturbation has been used, 
which however neglects the interband scattering. [T91I25] 
Using Eq. ([TS]) , the square of the matrix element for edge 
roughness reads 



|«fc'|C/|n,fc)|2 = 



AA2 



E'^ 



1 



2 \ ^ 2 



k h 



l,n un h,n 

"'F c _ 2 /2^2_£. 

K El 



(21) 



Regarding the diagonal contributions of the transition rate 
matrix, it is convenient to write the inverse of the relax- 
ation time for the nth subband as 
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For a large number of occupied subbands [M ^ 1), the 
second term in the bracket can be neglected since k^i,^ 
Ep, and the summation can be replaced by an integral. 
Ea. (P^ can then be written as 
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which shows a striking similarity for the corresponding 
results reported for semiconductor quantum wires. |15| . 
Furthermore, this relaxation time results in a conductivity 
in the limit of A/" ^ 1 of 
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where we have used the Gaussian integral in the evaluation 
of part in the y-direction. In the case of small correlation 
length A <C A^, the transition matrix element has the 
form 



Up to now, we have looked at the transport with edge 
roughness in the absence of magnetic fields B and now 
continue by including it and discussing its semiclassical 
effects. A weak magnetic field may homogenize the con- 
tributions of the occupied subbands to the overall conduc- 
tivity and result in a reduction of the magnetoconductiv- 
ity when the cyclotron radius ^ W . This effect has 
been verified numerically by recent quantum calculations 
in GNRs [14] but, to the best of our knowledge, has so 
far not yet been observed. Here the maximum reduction 
of the relaxation time t in magnetic fields can be roughly 
estimated by averaging over all occupied modes, resulting 
in 

^ « AA^vpEl. (25) 

The magnetoconductivity near the dip is then given by 
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which holds for 1 . 
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This results in a conductivity given by 
" - h ^9/2 AA2^/^^^ 

n,n 
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Numerical results and discussion. In Fig. [T]we 
show the conductivity for (5-typc impurities according to 
Eq. (|15p as a function of Fermi energy. For different de- 
grees of disorder, the conductivity shows very similar fea- 
tures while the amplitude of the conductivity depends on 
the disorder parameters. Prominent quantum oscillations 
at zero temperature are observed, i.e., the conductivity 
drops rapidly as a new scattering channel is opened and 
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Fig. 1; (Colour online) The numerical conductivity of 5- type 
impurities in armchair GNRs with W — ISOnm plotted versus 
the Fermi energy for two different temperatures T. We have 
chosen the parameters as 7 = 0.9 and rii = 0.5, where 7 = 
j{hvF)'^ / Ef and = ni\%. 

increases again until the Fermi energy hits the next sub- 
band at larger energies. The oscillations are smeared by 
finite temperature to some extent. In the whole range of 
Fermi energies studied, the average conductivity remains 
independent of the carrier concentration, which is consis- 
tent with the two-dimensional case. [2l[3] 

Fig. [2] shows the conductivity for edge roughness as 
a function of Fermi energy calculated from Eq. ([21]) . 
The parameter values chosen for A and A correspond 
to short-range defects, for instance, a few atoms miss- 
ing at the GNR edges, as widely assumed in simulations 
of edge disorder [8j[9l[T4j[26]. The correlation length en- 
sures that A ^ Xp over the whole range of Fermi ener- 
gies. The Boltzmann conductivity at nonzero tempera- 
ture (indicated by the solid lines) shows suppressed quan- 
tum fluctuations in comparison with the zero-temperature 
cases (indicated by the dashed line). The overall conduc- 
tivity decreases as the Fermi energy increases. Since the 
correlation length A and edge position fluctuation ampli- 
tude A increase relative to the Fermi wavelength as Ep 
is increased, this behavior is similar to that one found in 
conventional quantum wires |15[|17| . In the case of large 
number of subbands, Af ^ 1, the results from Eq. ([24]) 
(indicated by solid lines) exhibit the same overall trends 
and agree well with the exact ones except the absence of 
the quantum oscillations. 

The Fermi energies in Fig. [T] and [5] correspond to num- 
bers of subbands between 10 and 30. For smaller Fermi 
energies, i.e. a few occupied modes, conductivity shows 
more prominent quantum fluctuations and may deviate 
considerably from the asymptotic expressions. Moreover, 
it should be noted that our calculations based on the 
Boltzmann approach is valid for the case where the in- 
terband scattering is rather strong. This is the case when 
the mean free path is considerably shorter than the length 
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Fig. 2: (Colour online) The Boltzmann conductivity as a func- 
tion of the Fermi energy for different edge roughness in AGNRs. 
The dashed and solid lines correspond to zero temperature and 
a temperature of T = lOK, respectively. The smooth solid lines 
are calculated from Ea. (|24l) in the limit of A/" ^ 1. 

of the graphene ribbon. 

Fig. [3] shows the Boltzmann conductivity as a function 
of the GNR width for different edge roughness parameters. 
Here, only results for T > are presented. The overall 
conductivity for two roughness levels exhibits a parabolic 
dependence on the width, superimposed by quantum oscil- 
lations. This quadratic behavior may be seen more clearly 
from the analytical expression Eq. (P^ . as illustrated by 
the dashed lines. 

In the following, we give a rough estimate of the GNR 
conductivity in magnetic fields with amplitudes close to 
the position of the ERID, i.e., Tc ~ W. A more exact cal- 
culation would have to rely upon a calculation of the wave 
functions in magnetic fields, which can be obtained by 
solving the eigenequation of the Dirac Hamiltonian with 
magnetic fields included. [27l[28] This, however, should 
have only a marginal effect and we limit ourselves to the 
qualitative properties of the system close to the ERID. 

In Fig. [4l the conductivity at T = lOK is shown in 
the vicinity of the ERID for different Fermi energies. The 
parameters for edge roughness are fixed to A = 0.3nm and 
A = 6nm. For comparison, the corresponding zero-field 
conductivity with Ep = 200meV is plotted as well. The 
conductivity around the ERID increases linearly with the 
GNR width, in contrast to the parabolic dependence in the 
absence of a magnetic field. This linear relationship can 
be easily identified from Eq. ([25)1 and is also illustrated 
by the dashed lines in Fig. 2] As a consequence of the 
distinctly different dependencies of a on W, the ERID 
can be expected to be more pronounced in wider GNRs. 
This feature is in qualitative agreement with our previous 
quantum simulations. [14j . 

We conclude this analysis by commenting on the ob- 
servability of the ERID in realistic GNRs. First of all, the 
length of the GNR is irrelevant in the present treatment 
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Fig. 3: (Colour online) The width dependence of the conductiv- 
ity for different edge roughness in AGNRs with Ef = 200meV 
and temperature T — lOK. The solid lines show the exact con- 
ductivity from Eq. (|11|I by the T matrix inversion, and dashed 
lines correspond to the limit A/" 2> 1 from Ea. (|24|) . 

since diffusive transport has been assumed. Second, we 
have restricted ourselves to the case of rather small corre- 
lation lengths for the edge roughness. It is self-evident that 
a large correlation length suppresses the ERID in view 
the reduced diffusiveness of the scattering at the edges. 
Furthermore, the bulk disorder must remain at a suffi- 
ciently low level as indicated by the quantum simulations 
before, such that it does not mask completely the edge 
roughness scattering. We moreover expect that qualita- 
tively, the ERID does not depend much in the the type 
of edges, even though numerical simulations suggest that 
zigzag GNRs are more robust with respect to edge disor- 
der. [13] Similar analytical expressions for zigzag GNRs 
are possible in principle but more complicated due to the 
presence of surface states and the interdependence of the 
transverse and longitudinal wave vectors. 

In summary, we have studied the transport properties 
of AGNRs with short-range impurities and edge rough- 
ness within the framework given by the Boltzmann equa- 
tion. An edge-roughness-induced magnetoconductivity 
minimum suggested by the recent quantum calculations 
is confirmed by the Boltzmann results and should be- 
come observable experimentally if the correlation length 
of the edge roughness is not much larger than the Fermi 
wavelength and the bulk disorder is sufficiently low. It 
has been shown that the ERID induced by the magnetic- 
field-enhanced diffusive scattering at rough edges shows a 
behavior very similar to that one found in conventional 
semiconductor quantum wires, despite the fundamentally 
different energy dispersion. 

^ ^ ^ 
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Fig. 4: (Colour online) The magnetoconductivity around ERID 
as a function of the width for different Fermi energies at finite 
temperature T = lOK. (Note the logarithmic scale.) The 
zero-field conductivity with Ef = 200meV is also plotted at 
the top. The solid and dashed lines correspond to the results 
from Eg. lpTj) and EQ. (f26)) for A/" 2> 1, respectively, rc is the 
cyclotron radius. 
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